The combined Lagrangian advection method

We present and test a new hybrid numerical method for simulating layerwise-two-dimensional geophysical flows. The method radically extends the original Contour-Advective Semi-Lagrangian (CASL) algorithm by combining three computational elements for the advection of general tracers (e.g. potential vorticity, water vapor, etc.): (1) a pseudospectral method for large scales, (2) Lagrangian contours for intermediate to small scales, and (3) Lagrangian particles for the representation of general forcing and dissipation. The pseudo-spectral method is both efficient and highly accurate at large scales, while contour advection is efficient and accurate at small scales, allowing one to simulate extremely finescale structure well below the basic grid scale used to represent the velocity field. The particles allow one to efficiently incorporate general forcing and dissipation

Vortical control of forced two-dimensional turbulence

A new numerical technique for the simulation of forced two-dimensional turbulence (Dritschel and Fontane, 2010) is used to examine the validity of Kraichnan-Batchelor scaling laws at higher Reynolds number than previously accessible with classical pseudo-spectral methods,making use of large simulation ensembles to allow a detailed consideration of the inverse cascade in a quasi-steady state. Our results support the recent finding of Scott (2007), namely that when a direct enstrophy cascading range is well-represented numerically, a steeper energy spectrum proportional to k^(−2) is obtained in place of the classical k^(−5/3) prediction. It is further shown that this steep spectrum is associated with a faster growth of energy at large scales, scaling like t^(−1) rather than Kraichnan’s prediction of t^(−3/2). The deviation from Kraichnan’s theory is related to the emergence of a population of vortices that dominate the distribution of energy across scales, and whose number density and vorticity distribution with respect to vortex area are related to the shape of the enstrophy spectrum. An analytical model is proposed which closely matches the numerical spectra between the large scales and the forcing scale

Impeded inverse energy transfer in the Charney--Hasegawa--Mima model of quasi-geostrophic flows

The behaviour of turbulent flows within the single-layer quasi-geostrophic (Charney--Hasegawa--Mima) model is shown to be strongly dependent on the Rossby deformation wavenumber $\lambda$ (or free-surface elasticity). Herein, we derive a bound on the inverse energy transfer, specifically on the growth rate \d\ell/\dt of the characteristic length scale $\ell$ representing the energy centroid. It is found that \d\ell/\dt\le2\norm q_\infty/(\ell_s\lambda^2), where \norm q_\infty is the supremum of the potential vorticity and $\ell_s$ represents the potential enstrophy centroid of the reservoir, both invariant. This result implies that in the potential energy dominated regime ($\ell\ge\ell_s\gg\lambda^{-1}$), the inverse energy transfer is strongly impeded, in the sense that under the usual time scale no significant transfer of energy to larger scales occurs. The physical implication is that the elasticity of the free surface impedes turbulent energy transfer in wavenumber space, effectively rendering large-scale vortices long-lived and inactive. Results from numerical simulations of forced-dissipative turbulence confirm this prediction.Comment: 8 pages, 2 figures, accepted for publication in JF

Effective degrees of nonlinearity in a family of generalized models of two-dimensional turbulence

We study the small-scale behavior of generalized two-dimensional turbulence governed by a family of model equations, in which the active scalar $\theta=(-\Delta)^{\alpha/2}\psi$ is advected by the incompressible flow $\u=(-\psi_y,\psi_x)$. The dynamics of this family are characterized by the material conservation of $\theta$, whose variance is preferentially transferred to high wave numbers. As this transfer proceeds to ever-smaller scales, the gradient $\nabla\theta$ grows without bound. This growth is due to the stretching term $(\nabla\theta\cdot\nabla)\u$ whose ``effective degree of nonlinearity'' differs from one member of the family to another. This degree depends on the relation between the advecting flow $\u$ and the active scalar $\theta$ and is wide ranging, from approximately linear to highly superlinear. Linear dynamics are realized when $\nabla\u$ is a quantity of no smaller scales than $\theta$, so that it is insensitive to the direct transfer of the variance of $\theta$, which is nearly passively advected. This case corresponds to $\alpha\ge2$, for which the growth of $\nabla\theta$ is approximately exponential in time and non-accelerated. For $\alpha<2$, superlinear dynamics are realized as the direct transfer of entails a growth in \nabla\u, thereby enhancing the production of $\nabla\theta$. This superlinearity reaches the familiar quadratic nonlinearity of three-dimensional turbulence at $\alpha=1$ and surpasses that for $\alpha<1$. The usual vorticity equation ($\alpha=2$) is the border line, where \nabla\u and $\theta$ are of the same scale, separating the linear and nonlinear regimes of the small-scale dynamics. We discuss these regimes in detail, with an emphasis on the locality of the direct transfer.Comment: 6 journal pages, to appear in Physical Review

Dynamic Potential Vorticity Initialization and the Diagnosis of Mesoscale Motion

13 pages, 12 figures, 1 tableA new method for diagnosing the balanced three-dimensional velocity from a given density field in mesoscale oceanic flows is described. The method is referred to as dynamic potential vorticity initialization (PVI) and is based on the idea of letting the inertia-gravity waves produced by the initially imbalanced mass density and velocity fields develop and evolve in time while the balanced components of these fields adjust during the diagnostic period to a prescribed initial potential vorticity (PV) field. Technically this is achieved first by calculating the prescribed PV field from given density and geostrophic velocity fields; then the PV anomaly is multiplied by a simple time-dependent ramp function, initially zero but tending to unity over the diagnostic period. In this way, the PV anomaly builds up to the prescribed anomaly. During this time, the full three-dimensional primitive equations-except for the PV equation-are integrated for several inertial periods. At the end of the diagnostic period the density and velocity fields are found to adjust to the prescribed PV field and the approximate balanced vortical motion is obtained. This adjustment involves the generation and propagation of fast, small-amplitude inertia-gravity waves, which appear to have negligible impact on the final near-balanced motion. Several practical applications of this method are illustrated. The highly nonlinear, complex breakup of baroclinically unstable currents into eddies, fronts, and filamentary structures is examined. The capability of the method to generate the balanced three-dimensional motion is measured by analyzing the ageostrophic horizontal and vertical velocity-the latter is the velocity component most sensitive to initialization, and one for which a quasigeostrophic diagnostic solution is available for comparison purposes. The authors find that the diagnosed fields are closer to the actual fields than are either the geostrophic or the quasigeostrophic approximations. Dynamic PV initialization thus appears to be a promising way of improving the diagnosis of balanced mesoscale motions. © 2004 American Meteorological SocietySupport for this research has come from the Spanish Ministerio de Ciencia y Tecnología Grant REN2002-01343 and the program Ramón y Cajal 2001, as well as the U.K. National Environment Research Council (Grant GR3/11899)Peer Reviewe

Balanced solutions for an ellipsoidal vortex in a rotating stratified flow

Support for this research has come from the UK Engineering and Physical Sciences Research Council (grant number EP/H001794/1).We consider the motion of a single ellipsoidal vortex with uniform potential vorticity in a rotating stratified fluid at finite Rossby number . Building on previous solutions obtained under the quasi-geostrophic approximation (at first order in ), we obtain analytical solutions for the balanced part of the flow at . These solutions capture important ageostrophic effects giving rise to an asymmetry in the evolution of cyclonic and anticyclonic vortices. Previous work has shown that, if the velocity field induced by an ellipsoidal vortex only depends linearly on spatial coordinates inside the vortex, i.e. , then the dynamics reduces markedly to a simple matrix equation. The instantaneous vortex shape and orientation are encapsulated in a symmetric matrix , which is acted upon by the flow matrix to provide the vortex evolution. Under the quasi-geostrophic approximation, the flow matrix is determined by inverting the potential vorticity to obtain the streamfunction via Poisson's equation, which has a known analytical solution depending on elliptic integrals. Here we show that higher-order balanced solutions, up to second order in the Rossby number, can also be calculated analytically. However, in this case there is a vector potential that requires the solution of three Poisson equations for each of its components. The source terms for these equations are independent of spatial coordinates within the ellipsoid, depending only on the elliptic integrals solved at the leading, quasi-geostrophic order. Unlike the quasi-geostrophic case, these source terms do not in general vanish outside the ellipsoid and have an inordinately complicated dependence on spatial coordinates. In the special case of an ellipsoid whose axes are aligned with the coordinate axes, we are able to derive these source terms and obtain the full analytical solution to the three Poisson equations. However, if one considers the homogeneous case, whereby the outer source terms are neglected, one can obtain an approximate solution having a compact matrix form analogous to the leading-order quasi-geostrophic case. This approximate solution proves to be highly accurate for the general case of an arbitrarily oriented ellipsoid, as verified through comparisons of the solutions with solutions obtained from numerical simulations of an ellipsoid using an accurate nonlinear balance model, even at moderate Rossby numbers.PostprintPeer reviewe

The Roll-Up of Vorticity Strips on the Surface of a Sphere

We derive the conditions for the stability of strips or filaments of vorticity on the surface of a sphere. We find that the spherical results are surprisingly different from the planar ones, owing to the nature of the spherical geometry. Strips of vorticity on the surface of a sphere show a greater tendency to roll-up into vortices than do strips on a planar surface. The results are obtained by performing a linear stability analysis of the simplest, piecewise-constant vorticity configuration, namely a zonal band of uniform vorticity located in equilibrium between two latitudes. The presence of polar vortices is also considered, this having the effect of introducing adverse shear, a known stabilizing mechanism for planar flows. In several representative examples, the fully developed stages of the instabilities are illustrated by direct numerical simulation. The implication for planetary atmospheres is that barotropic flows on the sphere have a more pronounced tendency to produce small, long-lived vortices, especially in equatorial and mid-latitude regions, than was previously anticipated from the theoretical results for planar flows. Essentially, the curvature of the sphere's surface weakens the interaction between different parts of the flow, enabling these parts to behave in relative isolation

Wave and Vortex Dynamics on the Surface of a Sphere

Motivated by the observed potential vorticity structure of the stratospheric polar vortex, we study the dynamics of linear and nonlinear waves on a zonal vorticity interface in a two-dimensional barotropic flow on the surface of a sphere (interfacial Rossby waves). After reviewing the linear problem, we determine, with the help of an iterative scheme, the shapes of steadily propagating nonlinear waves; a stability analysis reveals that they are (nonlinearly) stable up to very large amplitude. We also consider multi-vortex equilibria on a sphere: we extend the results of Thompson (1883) and show that a (latitudinal) ring of point vortices is more unstable on the sphere than in the plane; notably, no more than three point vortices on the equator can be stable. We also determine the shapes of finite-area multi-vortex equilibria, and reveal additional modes of instability feeding off shape deformations which ultimately result in the complex merger of some or all of the vortices. We discuss two specific applications to geophysical flows: for conditions similar to those of the wintertime terrestrial stratosphere, we show that perturbations to a polar vortex with azimuthal wavenumber 3 are close to being stationary, and hence are likely to be resonant with the tropospheric wave forcing; this is often observed in high-resolution numerical simulations as well as in the ozone data. Secondly, we show that the linear dispersion relation for interfacial Rossby waves yields a good fit to the phase velocity of the waves observed on Saturn’s ‘ribbon’

Wave and Vortex Dynamics on the Surface of a Sphere

Motivated by the observed potential vorticity structure of the stratospheric polar vortex, we study the dynamics of linear and nonlinear waves on a zonal vorticity interface in a two-dimensional barotropic flow on the surface of a sphere (interfacial Rossby waves). After reviewing the linear problem, we determine, with the help of an iterative scheme, the shapes of steadily propagating nonlinear waves; a stability analysis reveals that they are (nonlinearly) stable up to very large amplitude. We also consider multi-vortex equilibria on a sphere: we extend the results of Thompson (1883) and show that a (latitudinal) ring of point vortices is more unstable on the sphere than in the plane; notably, no more than three point vortices on the equator can be stable. We also determine the shapes of finite-area multi-vortex equilibria, and reveal additional modes of instability feeding off shape deformations which ultimately result in the complex merger of some or all of the vortices. We discuss two specific applications to geophysical flows: for conditions similar to those of the wintertime terrestrial stratosphere, we show that perturbations to a polar vortex with azimuthal wavenumber 3 are close to being stationary, and hence are likely to be resonant with the tropospheric wave forcing; this is often observed in high-resolution numerical simulations as well as in the ozone data. Secondly, we show that the linear dispersion relation for interfacial Rossby waves yields a good fit to the phase velocity of the waves observed on Saturn’s ‘ribbon’

Toward a PV-based algorithm for the dynamical core of hydrostatic global models

The diabatic contour-advective semi-Lagrangian (DCASL) algorithms previously constructed for the shallow-water and multilayer Boussinesq primitive equations are extended to multilayer non-Boussinesq equations on the sphere using a hybrid terrain-following-isentropic (sigma-) vertical coordinate. It is shown that the DCASL algorithms face challenges beyond more conventional algorithms in that various types of damping, filtering, and regularization are required for computational stability, and the nonlinearity of the hydrostatic equation in the sigma- coordinate causes convergence problems with setting up a semi-implicit time-stepping scheme to reduce computational cost. The prognostic variables are an approximation to the Rossby-Ertel potential vorticity Q, a scaled pressure thickness, the horizontal divergence, and the surface potential temperature. Results from the DCASL algorithm in two formulations of the sigma- coordinate, differing only in the rate at which the vertical coordinate tends to with increasing height, are assessed using the baroclinic instability test case introduced by Jablonowski and Williamson in 2006. The assessment is based on comparisons with available reference solutions as well as results from two other algorithms derived from the DCASL algorithm: one with a semi-Lagrangian solution for Q and another with an Eulerian grid-based solution procedure with relative vorticity replacing Q as the prognostic variable. It is shown that at intermediate resolutions, results comparable to the reference solutions can be obtained.Publisher PDFPeer reviewe